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Abstract 

We analyze the prospects for discovering supersymmetry at the Fermilab Tevatron and 
CERN LEP colliders in the scenario that the lightest supersymmetric particle is a grav- 
itino of mass <, 1 keV. We consider in particular the case that the lightest neutralino has 
a nearly 100% branching fraction into gravitino + photon within the detector. This im- 
plies that supersymmetric events should contain both missing (transverse) energy and two 
energetic photons. Therefore one can search for supersymmetry simply through inclusive 
production of superpartners. We consider the exclusion and reach capabilities of the Teva- 
tron in exploring the supersymmetric parameter space, and study the efficiencies which can 
be achieved in this search. We also consider the discovery reach and backgrounds at LEP 
with yfi = 160, 175, and 190 GeV. 
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1. Introduction 



One of the intriguing theoretical aspects of super symmetry is that if it is realized as a local 
symmetry, it necessarily and automatically incorporates gravity. This connection is of no con- 
sequence in most studies of super symmetric phenomenology at colliders, because of the familiar 
negligibility of gravitational interactions. However, this need not be so if the gravitino (the 
spin 3/2 partner of the graviton) is very light. The gravitino (G) obtains its mass by absorbing 
the spin 1/2 would-be goldstino associated with the spontaneous breaking of super symmetry. 
In the high-energy limit, the interactions of the ±1/2 helicity components of the gravitino are 
the same as those of the goldstino it has absorbed. As emphasized originally by Fayet (l|, 
these interactions are proportional to 1/mg in the txIq — ► limit and are therefore potentially 
important even for processes at ordinary energies. 

However, the strength of gravitino interactions (or equivalently l/m^) certainly cannot 
be arbitrarily large. The gravitino mass is related to the scale of spontaneous supersymmetry 
breaking by 

-G = ^N = 5.9xW-^^) 2 eV. (1) 
G V3M V 500 GeV / 

Here M = (87rGNcwton) -1 ^ 2 = 2.4 x 10 18 GeV, and Ag^gy is the coupling of the would-be 
goldstino to the divergence of the supercurrent. Now, the scale Asusy should at least exceed 
the mass of the heaviest of the superpartners of the Standard Model (SM) particles, typically 
a gluino (g) or squark (q), and probably greatly so. [One might expect a significant hierarchy 
between Asusy and the electroweak scale, in order that negative radiative corrections to the 
Higgs scalar (mass) 2 can be effective in driving electroweak symmetry breaking.] Thus if one 
takes e.g. a bound Asusy > 500 GeV, equation ([j]) becomes a lower bound on the gravitino mass 
of roug hly 6 x 10~ 5 eV. In any case a given mass spectrum for the sparticles always implies a 
lower bound on m^. This is equivalent to a bound on the strength of the gravitino's interactions 
with the SM particles and their superpartners. This type of bound is quite conservative, and is 
certainly not expected to be saturated in particular models [||, ||] of super symmetry-breaking at 
low energies. For example, recently proposed models Q of dynamical supersymmetry breaking 
communicated to the visible sector by gauge interactions evidently favor tUq <; 1 eV which 
automatically avoids a dangerous R-axion [Q], although other ways of doing this might be 
possible. 

One can also attempt to obtain a lower bound on the gravitino mass by examining the 
requirements of partial- wave unitarity for, e.g., the scattering of two gluons into two gravitinos 
That process has contributions at tree-level from t- and u-channel exchanges of the gluino. 
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As shown in ref. || , tree-level partial wave unitarity is violated in this process when exceeds 

E cvit = V288irMm d /mg . (2) 



Now, one way to interpret this result is that the gravitino interactions, being gravitational, 
should not become strong below, say, the scale M; then (|2|) would become a lower bound on the 
gravitino mass of > m g /30. If this were true, the gravitino would always interact far too 
weakly to play any role in collider experiments. However, it seems preferable to interpret the 
critical energy indicated by (||) as the maximum value of a scale of unknown new physics 
at which m g becomes effective. In that case, one finds only that > A g m g /30M . If A g is 
smaller than the ultimate scale of supersymmetry breaking Asusy> this constraint is vacuous 
when compared with (jl]). 

On the other hand, cosmological constraints || seem to place an upper bound on tUq 
of about 10 4 eV, at least in the absence of late inflation. There is then still a window of 
perhaps 9 orders of magnitude for the mass of a light gravitino. In particular classes of models, 
this window can be much smaller. Throughout this window, m,Q is clearly insignificant for 
collider kinematics, and so can be taken to simply parameterize the strength of the gravitino's 
interactions. 

Most collider phenomenology studies performed up to now assume that the lightest super- 
symmetric particle (LSP) is a neutralino (mixture of neutral higgsinos and SU(2)l x U(1)y 
gauginos). If, as is most often assumed, i?-parity is exactly conserved, then supersymmetric 
particles will always be produced in pairs, and the LSP is absolutely stable. In this "neutralino 
LSP scenario", every supersymmetric event will feature two LSPs leaving the detector. There- 
fore the signals for supersymmetry always involve missing energy, often together with lepton 
and/or multi-jet signatures corresponding to particular decay chains of the superpartners pro- 
duced 0. In the "gravitino LSP scenario", however, the signatures should be quite different 
if the decays of superpartners into gravitinos occur within the detector a significant fraction of 
the time. 

For example, in the most obvious case that the next-to-lightest supersymmetric particle 
(NLSP) is a neutralino (iV"i) with a non-zero photino component, one has the interesting decay 

§§ 

Ni - 7 G . (3) 

In the rest frame of the decaying JVi, the photon takes an energy equal to m^j2 and is 
produced isotropically. The photons produced in supersymmetric events should therefore often 
be energetic enough to pass cuts designed to reduce SM backgrounds. The gravitino still carries 
away a significant amount of missing energy. Thus, supersymmetric signals in the gravitino LSP 
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scenario should be similar to those in the neutralino LSP scenario but with two (one) energetic, 
often isolated, photons if both (one) of the Ni decays occur within the detector. Several recent 
papers have presented interesting studies of this type of signal at the Next Linear Collider ||, 
ey colliders (l(|, and the Tevatron [11]. 

Since the presence of additional energetic photons in supersymmetric events would provide 
a welcome and powerful discriminant against SM backgrounds, it is possible to consider super- 
symmetry discovery signals based on inclusive production of all superpartners. The signatures 
in the gravitino LSP case are + X (or possibly + X), where X is an arbitrary col- 

lection of leptons+jets (including X = nothing which can occur e.g. in the cases of N\N\ or 
vv production). The purpose of the present paper is to study this strategy at the Fermilab 
Tevatron and CERN LEP colliders in the Minimal Supersymmetric Standard Model (MSSM) 
with a light gravitino. We define the MSSM to be the minimal supersymmetric extension of 
the SM with i?-parity conserved. Squarks (other than the top squarks) are assumed to be 
very nearly degenerate in mass, as suggested both by theory and the absence of flavor-changing 
neutral currents in experiment. Sleptons with the same electroweak quantum numbers are also 
assumed to be degenerate in mass. We will also often, but not always, make use of the "gaugino 
mass unification" assumption for running gaugino mass parameters: 

M 2 = r + \n M l = — M 3 (4) 

5 tan z U\y a 3 

which arises both in gauge-mediated and gravity-mediated supersymmetry breaking.^] 



This study is motivated in part by the observation at CDF of a single ee-y-yJ^T event [12] 



that does not seem to have a SM interpretation. It has already been suggested that this event 



can be explained by supersymmetry, either in the gravitino LSP scenario 13, 14 1 considered 



here or in a scenario with a higgsino-like neutralino LSP [14]. In [14|, we found that if this event 
is due to selectron pair production followed by the decays e — > eN\ and N\ — * jG, then the 
kinematic requirements of the event place rough bounds of 80 GeV < mg and 38 GeV < < 
100 GeV. However, this event can also be ascribed to pair production of charginos, as we shall 
remark below. In any case, we will maintain a more general point of view in most of the present 
work, rather than restrict our attention to the parameter space suggested by that one event. 

This paper is organized as follows. In section 2, we will consider the decays of supersym- 
metric particles into 2- and 3-body final states which include a gravitino. The 2-body decays 
N\ —* ZG and N\ — > h°G (which might compete with N\ — ► 7G if the photino component 
of Ni is very small) turn out to be subject to a very strong kinematic suppression. We also 
discuss expectations for the decay lengths of Ni, and note the existence of regions of parameter 

1 In particular, this "unification" relation can be well-motivated even in the absence of gauge-coupling unifi- 
cation at a very high energy scale; see e.g. fa]. 
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space where the decay N\ — > 7G is unduly suppressed. In section 3, we consider the limits 
on the reach of the Tevatron by studying the cross-sections for inclusive sparticle production. 
These rates are more interesting than in the case of the neutralino LSP scenario because of the 
relative ease with which arbitrary types of sparticle production can be detected using energetic 
photons. We propose a set of cuts designed to maximize the efficiency for detection of super- 
symmetry at the Tevatron via the signature + X, and study the efficiencies and lepton 
and jet multiplicities obtained using several sets of model parameters as test cases. In section 
4 we discuss the possibility of discovering supersymmetry at LEP with ^/s = 160, 175, and 190 
GeV, including an analysis of the relevant backgrounds. Section 5 contains some concluding 
remarks. 

2. Decays into gravitinos 

The decay Ni — > 7G will play the central role in the phenomenological discussions to follow. 
In sections 3 and 4 we will simply assume that this decay is the only important one involving the 
gravitino, and that it always occurs within the detector. However, it is interesting to consider 
other possible decays which might have an impact on collider phenomenology as well. If N\ 
is the NLSP and is fairly heavy, (but especially if it is beyond the reach of the Tevatron with 
the present integrated luminosity), the decays iVi — > h°G and Ni — > ZG can at least naively 
be important. Also it is interesting to consider the possibility that a sparticle other than N\ 
is the NLSP. Finally, if txIq is very small (<C 1 eV) one can even entertain the possibility that 
superpartners other than the NLSP can decay directly into final states containing a gravitino. 
In this section we present general formulas for the decay widths of supersymmetric particles 
into final states involving gravitinos. 

The relevant interactions of the gravitino are given by jl], |l5| 



contained in the field strength F^,, and (</>, ijj) are the scalar and fermionic components of 
the chiral supermultiplets. The full gravitino field can be well-approximated by its spin 1/2 
goldstino component when it appears as an external state in processes at energy scales relevant 
for collider studies: 



in the MSSM to final states including gravitinos. Let us first consider the decays of neutralinos. 





(6) 



In this limit it is not difficult to use eq. (||) to calculate decay rates of supersymmetric particles 
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Using the relation between the mass eigenstates and the gauge eigenstates of the one finds 
the decay widths of Ni into neutral gauge bosons to be as follows: 

771 ~ 



T(Ni -» ZG) = 2KiZ T + ^L ™\ L _ ™z_ | (8) 
v y 96vr M 2 m^, I m% 1 



1 



where 



Kj 7 = |iVii cos 9w + ^2 sin #vH (9) 
Kiz T = \Na sin 9w - N i2 cos 8w\ 2 (10) 
= \N i3 cos f3 - A^ 4 sin/3| 2 (11) 

measure the contents in iV, of photino, zino, and the higgsino partner of the neutral would-be 
Nambu-Goldstone boson, respectively. (Here and in the following we use the notations of |l(| 
for the parameters and mixing matrices of neutralinos and Higgs scalar bosons in the MSSM. 
Thus Nij are the neutralino mixing matrices with (i, j) the (mass, gauge) eigenstate labels, and 
tan/3 is the ratio of Higgs vacuum expectation values.) In these and similar formulas below, 
the Planck-scale suppression m 2 ^ /M 2 is numerically counteracted by the hierarchy m 2 ^ Jm 2 ^, 
so that the decay width can be non-negligible. So for example we can write © in the suggestive 
form 

/ 711 " \ / Til ~ \ — ^ 

r(*.-»fl) = i.i2xur»G.v^( iSi ^) (^) . (12) 

When the 2-body decay Ni — > ZG is near threshold, the formula (||) is not reliable; also 
when the 2-body decay is not kinematically allowed, the decay can still proceed through an 
off-shell Z boson. In these situations one must use the 3-body decay formula. In the limit 
of massless SM fermions from the off-shell Z, the width of the neutralino from 3-body decays 
through a virtual Z boson are obtained by replacing 

(tti 2 \ 
1 j- — > lKiZ T h + KiZ L Io (13) 
m Nj 

in (||), where the kinematic factors are most compactly written as 

< i4> 

7T Jo (x — R) z + e z 

with R = nig/m 2 ^ and e = Tzmz/rn 2 ^ . In the case that rn^ — mz 3> Tz, one finds Iq rj I\ « 



(' 



l 



1 — m z /m ~ \ so that (||) is recovered. At threshold (m^_ ~ m^) one finds Iq ~ 4J\ ~ .0029, 
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rather than as suggested by the 2-body formula (||). For another point of reference, with 
= 70 GeV one finds the kinematic factors are approximately Iq ~ .0013 and I\ ~ .00015. 
In order for off-shell or near-threshold decays N% — > Z^G to compete with N\ —* jG, the 
photino component of iVi clearly would have to be very small. Even for = 150 GeV one 
finds J » h » 0.16. 

Next we consider the decays of a neutralino into a gravitino plus any of the neutral Higgs 
scalar boson mass eigenstates (p = (h°,H°,A°) of the MSSM. The 2-body decay widths are 
given by 



in 



5 / 2 



r(JVi -> <pG) = gjjjj^r y ~ ^f- 1 («) 

where the relevant higgsino contents are given by 

K ih o = | Nis sin q — A^j4Cosa| 2 (16) 
Kj#o = |A^j3 cosa + iVj4sina| 2 (17) 
«Uo = |iV i3 sin/3 + iV i4 cos^| 2 (18) 

Note that with the identifications (p = G° (the electroweak would-be Nambu-Goldstone boson), 
K iG° = \Ni3c0sP — iVj 4 sin/3| 2 and m G o = mz, one recovers the decay width into longitudinal 
Z's indicated in equation (||), in compliance with the equivalence theorem [17]. 



It is certainly not an outlandish possibility that rrijj > m h o, so that the two body decay 
N\ — > h°G can compete with iVi — > 7G. However, this decay is also crippled by the kinematic 
suppression indicated in ( p"5| ) unless "m^ is significantly larger than m h o. The 3-body decay 
widths for N{ — > ip*G with ip* decaying into pairs of SM fermions (treated as massless for 



kinematic purposes) are given by replacing in (15) 

,2 



l--f (19) 



as given by (14) with R = m 2 /m 2 ~ and e = T^m^/m 2 - . Since in the MSSM the width of h 



is only a few MeV, such three-body and near-threshold decays are generally negligible. 

The decay widths of sleptons and heavy squarks are also easily found. For a sfermion 
decaying into a massless SM fermion + gravitino, one finds the 2-body width 

r ^' e ) = «fe (20 > 

Cr 

One of the more intriguing possibilities is that the nearly degenerate right-handed sleptons 
e R , fl R , f R act effectively as co-NLSPs. In that case, the supersymmetry discovery signatures 



7 



would generally involve at least two energetic leptons + $t- If a sneutrino is the NLSP, then 
signatures should often be similar to those in the neutralino LSP scenario, since the decays 
v — * vG are invisible. It does not appear to be sensible to contemplate a left-handed charged 
slepton as the NLSP, because of the sum rule 



m~ = m~ + m 2 v \ cos 2/3| (21) 



for tan/? > 1. Squark (mass) 2 parameters receive large positive contributions proportional 
to a 2 and/or a s , so that it seems doubtful that a squark could be the NLSP. One possible 
exception is the lightest top squark mass eigenstate (ti) since m~ can receive large negative 
radiative corrections proportional to the top Yukawa coupling squared. However, if < m t , 
then t\ should be very long-lived if it is the NLSP, and in particular should always hadronize 
and escape the detector as a charge or 1 "mesino" (iiq) or as a charge 0, 1, or 2 "sbaryon" 
{t\qq') bound state. In any case, for the remainder of this paper, we will decline to consider 
the possibility that a sfermion could be the NLSP. 

The 2-body decay widths of charginos (Cj) into gravitino final states are given by formulas 
entirely analogous to (JsJ) and (|l5|): 



r [Ct -+w G)- — I i — i (-) 




r(C+ ^ H+G) = KiH+ *\ I (23) 

v 1 96?r M 2 m| I m 2 d ' 



KiW T = U\Vii\ 2 + \U il \ 2 ) (24) 



with 



2 

Kiw L = \V i2 \ 2 sin 2 (3+\U i2 \ 2 cos 2 (25) 
K iH + = \V i2 \ 2 cos 2 p + \U i2 \ 2 sin 2 p (26) 

The generalizations to off-shell decays are given by the obvious analogs of the above expressions 
for Ni decays. However, it should be noted that because of the form of the chargino and 
neutralino mass matrices, a chargino can only be the NLSP in a small and not particularly 
attractive region of parameter space. 

In general, if could be arbitrarily small compared to superpartner masses, then all 
decays of supersymmetric particles could proceed directly to the corresponding SM particle 
plus gravitino. However, as a practical matter for supersymmetric states accessible at Tevatron 
energies and taking into account a conservative lower bound on the gravitino mass as mentioned 
in the Introduction or stricter bounds in particular classes of models, it is easy to see that 
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the decay widths for non-NLSP sparticles listed above should be quite small and should be 
overwhelmed by the usual well-studied decays. Besides the decay of the NLSP, there is one other 
potential exception which seems worthy of mention. If the gravitino mass is near the lower end 
of the allowed window, it is possible that a heavy gluino can decay directly to gluon+gravitino 
through a 2-body decay, rather than following the usual cascade decay pattern through virtual 
squarks. Since the only other decays of the gluino are mediated by virtual squarks which can be 
quite heavy in models, it is conceivable that the direct decay to gravitino can dominate in the 
gravitino LSP scenario. (In contrast, decays of all other non-NLSP superpartners can proceed 
through virtual Ws or Zs or sparticles which are plausibly much lighter than squarks.) The 
relevant decay width is given by 

T (g _> g G) = ^L- = 1.1 x io-9 GeV ( m l V ( (27) 

vy y ' 487rM 2 m% V250 GeVy \l eV J V ; 

The competition between this decay and the usual cascade decays of gluinos has already been 
studied in |l8|, where (|?]) was found to be negligible unless m,Q < 10~ 2 eV, even if all 
squarks are as heavy as several TeV. For gravitinos lighter than 10 -2 eV, it was found that 
the direct decays ( p7|) can dominate over the more conventional decay chains through virtual 
squarks only if 

m o * 10 ~ 3eV (iiv) 2 - (28) 

This can only occur in the slightly problematic case that Astjsy does not greatly exceed rriq [cf. 
equation ([[])]. 

Given the considerations above, we will optimistically assume for the remainder of this 
paper that a neutralino is the NLSP, that the branching fraction for N\ — ► jG is 100%, and 
that all super symmetric decay chains terminate in this subdecay. Nevertheless, it is worth 
noting that this decay can be strongly suppressed due to a very small photino content of N\ in 
regions of parameter space with small making N\ long-lived on collider scales. Assuming 
the usual gaugino mass unification condition (|j) and restricting our attention to the parameter 
space not already excluded by LEP with tan/3 > 1.5, ki 7 can be less than 0.001 only if /i is 
negative and \fi\/M2 < 0.2. For « l7 < 0.01, it is required that \fj,\/M2 < 0.4 (0.2) for fi negative 
(positive). A milder but still quite significant suppression ki 7 < 0.1 can be obtained if \/j,\/M2 < 
0.5 (0.65) for \x negative (positive). Conversely, as long as > M 2 and > 50 GeV, one 
finds Ki 7 > 0.21 (0.13) for [i negative (positive). Formally, if k\~ were to vanish, the decay 
N\ — > jG could still proceed through one- loop graphs, but these amplitudes are very small |2(| 
and in the present context are only competitive for a fine-tuning of ki 7 which is finer than we 
are willing to contemplate here. 

Of course, a 100% branching fraction for Ni — > jG does not guarantee that the photons 
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can be detected, since the length scale associated with this decay might easily be comparable 
to the relevant physical size of the detector. The probability that each N\ with energy E in 
the lab frame will travel a distance < x before decaying is given by 

P{x) = 1 - e~ x/L (29) 

where from ([/]) the decay length is 

l _ i.76 x io- w -> ( £ >4, - D 1/2 (^kvY' (B) 2 cm • (30) 

Clearly L depends strongly on and tjIq and can either be larger than, or negligible compared 
to, the relevant physical dimension (~ 150 cm) of a CDF-type detector. Note that if L is larger 
than 150 cm, the efficiency for detecting one photon can greatly exceed that for detecting both. 
For example, taking L to be 15 meters, one finds that the probability for both (one) of the 
photons being emitted within 150 cm of the event vertex is roughly 0.01 (0.17). For L ~ 150 
cm, the probability of two (one) photons being emitted within 150 cm of the event vertex is 0.40 
(0.47). Since the SM backgrounds for events with one energetic photon greatly exceed those 
for events with two such photons, we will optimistically assume in the following discussion that 
for the processes of interest L < 150 cm, so that all supersymmetric events will lead to two 
potentially detectable photons. Taking ki 7 and (E 2 /m 2 ^ — l) 1 / 2 to be of order unity, this 
requires roughly rriQ < 250 eV for a 100 GeV N\. Larger decay lengths will decrease the 
efficiency of detection accordingly. 



3. Supersymmetry with a light gravitino at the Tevatron 



The presence of two energetic photons from supersymmetric events in the gravitino LSP 
scenario should dramatically increase the detectability over that found in the usual neutralino 
LSP scenario. In this section we will study the possibility for detecting inclusive + 
X signals at the Tevatron in the present data sample of about 100 pb _1 per detector. We 
concentrate on signals from chargino and neutralino, slepton, and light stop squark production 
for a range of models, and we comment on other potential signals. For this study, we assume 
that the decays N\ — > 7G occur within the detector 100% of the time. As a practical matter, 
we compute kinematics of events with the further assumption that these decays occur close 
to the event vertex. All event simulation is performed using the PYTHIA Monte Carlo with 
supersymmetric extensions Q|. Refs. || [TI], [ll[] also contain recent studies of gravitino LSP 
physics at colliders. 
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A. Chargino and Neutralino Production 

The production cross sections for CiCj, CiNj and NiNj at hadron colliders are functions 
of the gaugino-higgsino parameters (the U{\) gaugino mass M\, the SU(2) gaugino mass M2, 
the supersymmetric Higgsino mass parameter fj,, tan/3) and the squark masses. In the following 
we vary these parameters to find the range of expected signals. The gravitino LSP scenario has 
striking phenomenological implications for Cj and iVj production. Consider, for example, the 
process pp — > C^iV^. In the neutralino LSP scenario, it is well-known that this process can 
be detected with a good efficiency when the final state includes three leptons p2| . However, 
in the gravitino LSP scenario, all of the final states of this process (including r y^£ + £^£' + f!T, 
77^ + jjJpT-, 77-^ '£~ jj$T and 'y^jjjj^r) can provide useful signals. Likewise, the production 
of CfCi pairs can lead to observable signals "fj£ + £'~ $t, ll^jj$T, and yyjjjj^T- When 
the gaugino unification condition (|j) is satisfied, these two processes provide the bulk of the 
supersymmetric signal throughout much of parameter space, because of the relatively large 
couplings WC1N2 and ZC\Ci, 7C1C1. Also N\N\ production, which is undetectable at hadron 
colliders in the neutralino LSP scenario, leads to the signal 77-^T in the gravitino LSP scenario. 
Unfortunately, although this process is kinematically favored, it usually has a negligible cross 
section at hadron colliders because of a small ZN\N\ coupling and heavy squarks. 

The branching fractions for the various final states associated with chargino and neutralino 
production are quite model-dependent. For example, the jet or lepton multiplicity and kine- 
matics can be a strong function of squark and slepton masses. Since all of these final states 
from Ci and N production involve two energetic photons and missing transverse energy, we 
prefer to focus on the inclusive 77^7^ + X signal rather than details of jet or lepton multiplicity. 
If a number of events are found in the data sample which are not understood as coming from 
the SM, then such details could help disentangle the underlying theory. Below we show such 
distributions for a specific set of models. 

In Fig. we show the allowed range for the total inclusive production cross-section a(pp — > 
NiNj or CiCj or NiCj) at ^/s = 1.8 TeV as a function of the lightest neutralino mass rn^ . 
We have assumed that the gaugino mass unification assumption @ holds, so that the gaugino- 
higgsino sector is determined by only three parameters, one of which we choose to be the N\ 
mass. To generate this graph, we have varied the other parameters of the MSSM over the 
ranges 

250 GeV < rriq < 1000 GeV 
-1000 GeV < n < 1000 GeV 
1.5 < tan/3 < 55 
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N, mass (GeV) 

Figure 1: Total production cross sections for charginos and neutralinos (NiNj and NiCj and 
CiCj) at the Tevatron with ^/s = 1.8 TeV as a function of the lightest neutralino mass, assuming 
gaugino mass unification. The solid lines are the minimum and maximum allowed cross sections. 
The dashed line is a typical large \fi\ and heavy squark limit (fj, = 1000 GeV, mq = 1000 GeV, 
tan/3 = 1.5). 

The dashed line represents a typical large heavy squark limit, namely [i = 1000 GeV, 
nig = 1000 GeV, tan/3 = 1.5. If the gaugino mass unification condition (||) is not satisfied, then 
the total inclusive chargino- neutralino production cross-section (for the range of rrifj shown) 
can be essentially negligible; this is traceable directly to the kinematic suppression associated 
with very heavy charginos. In Fig. || we show the same cross-section, this time as a function 
of the lighter chargino mass. Again, Fig. § assumes (|j). However, we found that the minimum 
production cross-section in the case of general gaugino mass parameters is not significantly 
lower than that shown in Fig. ^, for a given C\ mass. This is important because it shows that 
the Tevatron can set model-independent exclusion limits on , if the efficiency for detection 
is reasonably bounded from below. The maximum production cross-section for general gaugino 
mass parameters not obeying (|j) can be several times larger than that shown in Fig. ^. As 
can be seen from Fig. [2|, the total number of chargino/neutralino pair production events at the 
Tevatron with 100 pb _1 of data is > 100 before cuts if a chargino mass is less than 100 GeV 
(the maximum LEP reach). 

To fully define a signal, we choose the following cuts for the two photons 71,72: 

• Et(ii), Ex{ r )2) > 12 or 25 GeV, where Et is the transverse energy. We define two 
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C, mass (GeV) 



Figure 2: As in Fig. |T|, but as a function of mg . We find that even when the gaugino mass 
unification assumption (0) is not made, the minimum cross section is never significantly less 
than shown here. 

different signals based on the minimal Et- For the lower threshold, we impose a 30% loss of 
efficiency when one or both photons have Et below 25 GeV, to simulate the approximate loss 
of triggering efficiency p3|| , 

• I ^(Ti ) I ? 1 77(72) | < 1) where 77 is the pseudorapidity. 

• £f < 4 GeV 0, where 

E™ = ]T E^-E T (n), 

j,R<0A 

and we sum the transverse energy from all particles (j) within a cone of size R = \J (Ar/) 2 + (A(p) 2 
centered on the photon candidate. (Photons from jets or bremsstrahlung tend not to be isolated 
from additional hadronic activity.) 

• $t > 30 GeV, where $t is determined by the sum of the visible energy in smeared jets, 
photons, and leptons. 

Standard Model physics backgrounds can arise from W {— > £1^)77, Z(—> vi?,T + T~) r yj, 
and QQ77 where Q = c,b, or t. We have not made a full matrix element simulation of these 
backgrounds. However, we estimate their magnitude by generating W^j, Zj, ti, bjX, and 
2 This number can be significantly reduced without affecting the signal. 
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cjX events with additional QED radiation in the leading-logarithm approximation [24|. Based 
on this analysis, we estimate (0.13, 0.11, < 10~ 3 ,< 10~ 2 ,< 0.1) events in 100 pb _1 from each 
of these sources when Et(j) > 12 GeV. For this set of cuts, the C77 background is the largest 
hadronic source, and it is well under control. Likewise, backgrounds from a jet faking an isolated 
photon can be estimated from these numbers using a simple scaling by a s /a em x Rj^-y X / ~ .1, 
where Rj^ ~ 10 -3 is the probability a jet fakes a photon and / ~ 10 accounts for the squared 
quark charge. The probability of two jets faking two photons is even further suppressed. Finally, 
backgrounds with a fake missing transverse energy are limited by excellent electromagnetic 
calorimetry. Essentially, these chosen cuts should yield a signal free from background, though 
the fix cut could be increased if necessary. 

To assess the sensitivity of the Tevatron to 77^+ X signals from C{ and Ni production, we 
have performed event level simulations of various light gravitino models using the parameters: 

M 2 = 100, 150, 200, 225 GeV; 
\x = ±125, ±215, ±300, ±600 GeV; 
tan/3 = 1.5, 1.7, 2.5, 3.0, 10.0; ^ > 

(mq,m lL ,m lR ) = (250,125,119) or (500,250,238) GeV. 

Here Mi is fixed by (Q). The sneutrino mass is fixed by the sum rule (|2l|). (When the result 
is less than the Aq mass, we take instead rriy = rrif[ 1 ±5 GeV and mj £ fixed by the sum 
rule.) While larger squark masses can easily be obtained in models, we find that to a good 
approximation, the dependence of signals on squark mass vanishes for squark masses above 
500 GeV at the Tevatron. We display the results in terms of the number of expected events 
in 100 pb -1 as a function of C\ mass in Figs. || and ||, for a minimum .Br (7) threshold of 12 
and 25 GeV respectively. The efficiency for detection of the signal is also displayed in Fig. [| 
as a function of and in Fig. | as a function of mg using the -&r(7) > 12 GeV cut. By 
comparing Figs. || and |], we conclude that most of the photons originating from models with 
m Ci — GeV which pass the 12 GeV Er cut will also pass the higher threshold. It is also 
clear that the lower £t(7) threshold substantially increases the signal for smaller despite 
the loss of triggering efficiency. 

These figures suggest that the non-observation of a signal excludes mg < 125 GeV when 
(Q) is assumed, which is well above the pair production threshold of any LEP upgrade. The 
same information, but plotted as a function of the N% mass, is shown in Fig. |?j. From this 
plot, we conclude that a lightest neutralino mass below 70 GeV is excludable in the same 
manner, when @) is assumed. We have not attempted a completely general study of efficiencies 
when the gaugino mass unification is not assumed. However, we do not see any reason to 
expect significantly lower efficiencies in the completely general case. In particular, small mass 
differences between charginos and neutralinos should have little effect on the efficiency (for 
fixed mrifj ) since the photon energies and J^t, which primarily determine the signal, depend on 
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Figure 3: The number of + X events expected in 100 pb _1 of data using the models 

defined by ([H]) and the cuts explained in the text with Et{j) > 12 GeV. The signal comes 
from the inclusive production cross sections for charginos and neutralinos {NiNj and NiCj and 
CiCj) at the Tevatron with y/s = 1.8 TeV as a function of the lightest chargino C\ mass, 
assuming gaugino mass unification. 
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Figure 4: As in Fig. |, but with E T (j) > 25 GeV. 
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Figure 5: The number of generated + X events passing cuts divided by the total for 

Et(j) > 12 GeV, for the models defined by (|3l|). The signal comes from the inclusive production 
cross sections for charginos and neutralinos (N{Nj and iVjCj and C\Cj) at the Tevatron with 
yfs = 1.8 TeV as a function of the lightest chargino C\ mass, assuming gaugino mass unification. 
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Figure 6: As in Figure ||, but as a function of tuq^. 
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Figure 7: The number of +X events expected in 100 pb _1 of data using the cuts explained 
in the text with £^(7) > 12 GeV. The signal comes from the inclusive production cross sections 
for charginos and neutralinos (NjNj and N^Cj and CiCj) at the Tevatron with y/s = 1.8 TeV 
as a function of the lightest neutralino N% mass, assuming gaugino mass unification, and using 
the models defined by (|3l]). 

the mass and boost of N\. Therefore by considering Figs. [5] and |6] and using the fact that the 
minimum production cross-section as a function of is bounded from below as in Fig. ||, we 
conclude that it should be possible to exclude < 100 GeV for > 50 GeV using the 
present 100 pb _1 of Tevatron data, even without assuming (Q). 

As mentioned previously, the lepton and jet multiplicities of such events can be large, 
although they can be sharply reduced from naive expectations because of limited detector 
acceptance, jet definition, and isolation criteria. This is particularly relevant when the mass 
splittings among charginos and neutralinos are relatively small. Jets (j) are defined using a 
standard clustering algorithm with R=0.5 and E 3 T > 15 GeV, |r/ J | < 2.5. The particle energies 
are smeared using typical CDF energy resolutions. Electrons and muons must have E^ > 20 
GeV and |r/ e '^| < 2.0, while being isolated from excess transverse energy. We illustrate typical 
jet and lepton multiplicities for four specific models in Fig. 0. The model parameters are: 



17 




2 4 

Multiplicity for P} > 25 GeV 



Figure 8: The lepton and jet multiplicities for the four models explained in the text. 

• Model 1: M 2 = 100 GeV, n = -216 GeV, tan/3 = 2.5, m$ = mj = 1000 GeV. One then 
finds mfj. = (53, 108,227,238) GeV, mg. = (108,240) GeV, «i 7 = 0.85. 

• Model 2: M 2 = 150 GeV, n = -125 GeV, tan/3 = 1.7, m$ = mj= 1000 GeV. One then 
finds mft. = (80, 119, 153, 179) GeV, m & . = (134, 182) GeV, «i 7 = 0.81. 

• Model 3: M 2 = 200 GeV, u = 600 GeV, tan/3 = 3.0, m~ p = 150 GeV, m = 250 
GeV, m ih = 260 GeV, m 5 = 700 GeV. One then finds m^. = (98,190,602,615) GeV and 
m di = (189, 614) GeV and « l7 = 0.72. 

• Model 4: M 2 = 225 GeV, n = 300 GeV, tan/3 = 1.5, = 105 GeV, m p = 115 
GeV, m ih = 125 GeV, m~ q = 300 GeV. One then finds = (101,183,301,355) GeV and 
m^. = (176, 350) GeV and «i 7 = 0.56. 

We note in passing that Model 4 may be of particular interest, since it has some general 
properties consistent with an alternative candidate for the CDF eejj^r event through CiCi 
production. In this model, the lepton multiplicity is peaked at 1 but there is also a substantial 
component with lepton multiplicity 2 (because of allowed 2-body decays of C\ and iV 2 to slepton 
+ N\). One expects about three chargino/neutralino events after cuts from this model in the 
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Fi gure 9: Total slepton production cross sections at the Tevatron with y/s — 1.8 TeV, for l R l R 
as a function of (solid line), and left-handed sleptons (pt> and l L v and £ L £ L ) as a function 
of rriy (dashed lines). In the latter case, the lower (upper) dashed line corresponds to tan (5 = 55 
(1.5). The cross sections shown are summed over slepton flavors, with slepton masses taken to 
be flavor- independent. 

current data taken at the Tevatron. We will remark further on the chargino pair production 
interpretation of the CDF event below. This model also has light sleptons, so it could have 
produced the event through selectron pair production, but the kinematics do not favor this 
interpretation because the leptons would be too soft. Fig. || shows how the relative multiplicities 
could help distinguish models if a signal is established. 

B. Sleptons 

In most theoretical models, scalar (mass) 2 parameters receive positive contributions pro- 
portional to af and/or aj, where a, [i = 1,2,3) are the gauge couplings felt by the scalar. 
Therefore one expects that sleptons with the same SU{2)l x U(1)y quantum numbers should be 
degenerate in mass, and should all be considerably lighter than squarks, with rrig < < 
the most plausible mass ordering. It is therefore interesting to consider slepton discovery sig- 
nals at the Tevatron; a corresponding study in the neutralino LSP scenario appears in pq| . In 
Fig. ^, we show total Tevatron cross-sections for i R i R production summed over three families 
(solid line) as a function of . The signal for e R e R production with the decay e R — ► eN\ 
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Figure 10: Expected £ R £ R signals (including all three lepton flavors) from Model 4 at the 
Tevatron with ^/s = 1.8 TeV, as a function of . The total 77 $r + X signal is shown by the 
filled circles, while the single lepton and dilepton components are denoted by crosses and open 
circles, respectively. 



is eej'j^T, providing a viable candidate for the single observed CDF event of this type [12] 



As is known from the analyses of 12, 13, 14|, such events do not seem to have a probable SM 
interpretation. In the same figure we show as a function of mp the total cross-section for DD, 
£ L i), and t L t L production, for tan/3 = 1.5 and 55 (dashed lines). Since the masses of u and £ L 
are related by the sum rule (^T|), the rates for £ L v and £ L £ L production decrease monotonically 
with larger tan/3 (for a fixed value of mp). The £ L v component of the signal is always the 
largest. The final states from £ L £ L and £ L v production will depend specifically on the slepton 
and Q, Ni masses, but can contain and up to three charged leptons. 

Rather than conduct an extensive survey of slepton signatures, we consider as a test case the 
chargino/neutralino sector of Model 4 of subsection A. For the fixed set of gaugino parameters 
of that model, we further vary the right-handed selectron and sneutrino masses over the ranges 
100 GeV < m^^rriy < 200 GeV. It should be noted that for much of this range, direct decays 
of slepton to lepton and N\ should dominate for the chosen model. Since the signals from right- 
handed sleptons and from left-handed sleptons have rather different characteristics, and because 
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Figure 11: Total expected i L i L and t h v and vv signals (including all three lepton flavors) from 
Model 4 at the Tevatron with ^fs = 1.8 TeV, as a function of rriy. The total 77$r + X signal 
is shown by the filled circles, while the single lepton and dilepton components are denoted by 
crosses and open circles, respectively. 



the masses of right- and left-handed sleptons are a priori unrelated, we choose to display the 
results separately. In Fig. [h] we show the cross-sections after cuts arising from right-handed 
slepton production. The inclusive 77-^T signal, with or without additional leptons, is denoted 
by filled circles. We note that the dilepton component of the signal (open circles) is greater than 
the single lepton component (crosses) for wig > 130 GeV. There is also a significant component 



with no leptons passing the cuts. In Fig. 11 we likewise show the total rate after cuts expected 
from £ L £ L and £ L v, and vv production, as a function of m„. For any given model, the expected 
number of events with zero or one lepton far exceeds the number with two leptons. This is 
partly because of the comparatively larger cross-section for vl L and vv production, but also 
because some leptons do not pass cuts. 



C. Light stop squarks 



In specific models, one stop squark mass eigenstate (ti) is often found to be much lighter 
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Figure 12: Cross section for pair-production of the lighter top squark mass eigenstate at the 
Tevatron with y/s = 1.8 TeV. 

than all of the other squarks, and can even be lighter than the top quark. If < 100 GeV, 
chargino/stop loops might help to explain |^] the excess of Rb in the LEP data. However, 
we have already seen that in the gravitino LSP scenario, charginos must be far too heavy for 
this to occur. Furthermore, a significant bound can be independently placed on a light stop 
squark mass in the gravitino LSP scenario given the integrated luminosity already obtained at 



the Tevatron. In Fig. |12j, we show the total t\t\ pair production cross-section as a function of t\ 



mass. In the gravitino LSP scenario, this process should lead to spectacular signals 77$r+jets. 

We consider two scenarios, based on the mass orderings < +rrib or > +mb- 
In the first case, each t\ cascades through two two-body decays to a C7G final state. In the 
second, t\ undergoes an additional three-body decay to reach a 6//7G final state, where / is 
a fermion. As a result, the photons produced in the second case tend to be softer. For the 
first case, t\t\ production leads to two additional charm jets in the final state, while 6-jets and 
additional leptons or jets are present for the second. As before, we ignore such particulars, which 
could substantiate a suspected signal, and concentrate on the same inclusive 77-^T + X signal. 
Based on the previous bounds on mjj and Tn^, we consider models with tUq + m& > > 
rrift > 70 GeV for the range of gaugino parameters discussed previously, and a smaller set of 
models with > uiq +rrib > 125 GeV. The results are illustrated in Fig. which shows the 
number of expected diphoton events in the present data sample using the previously defined 
cuts with Exd) >12 GeV. There is a substantially higher detectability of the signal when 
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Figure 13: Number of + X events in 100 pb 1 as a function of the lighter top squark 

mass eigenstate at the Tevatron with yfs =1.8 TeV. 

t\ — > C\b is not kinematically allowed. The absence of such events in the present data sample 
seems to exclude a t\ lighter than at least 140 GeV even in the case that > fng +mb, which 
is already far too heavy to have any effect on the interpretation of the LEP Rb measurement. 
Additionally, if + < mj, then the decay t — > t\Ni can occur, generating 77-^r events 
from ti production, but the bounds on and mfi preclude this. 

Note that the limit on the mass of the light stop in the gravitino LSP scenario is much 
stronger than for the neutralino LSP scenario in the case where t\ — > cN\. The latter case, 
which relies on the signal of two acollinear jets and is limited by the mass splitting between 
t\ and Ni which determines the jet spectrum f27|| . The main limitation of the gravitino 
LSP scenario is which sets both the scale of Et{^) and $t- 

D. Other Processes 

It is a common feature of known models that the gluino and squarks are quite heavy. 
However, the presence of energetic photons in the eventual decay products means that the 
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Figure 14: Contours of inclusive production cross sections for gluinos and squarks (gg and 
gq, qq) at the Tevatron with y/s = 2 TeV. All squark flavors are taken to be degenerate for 
simplicity. 

detection efficiency is likely to be higher in the gravitino LSP scenario than in the neutralino 
LSP scenario. Therefore it is again interesting to get an idea of the upper limit on the potential 
reach of the Tevatron collider by considering the total inclusive production of gluinos and 
squarks (gg, gq, qq). It must be mentioned that there are at least two factors which might 
adversely affect the efficiency here somewhat, including larger boosts leading to longer decay 
lengths for N\ — > 7G [cf. eq. (|30|)1, and losses from photon isolation requirements with the 



higher jet multiplicity. In Fig. |14j, we show contours of this total production cross section, in 
the (mg,rriq) plane. (For simplicity, we have assumed degeneracy of all squark flavors.) At least 
in the case when gaugino mass unification is assumed, it appears to be doubtful that gluino 
pair-production processes can be useful, since the non-observation of events from chargino and 
neutralino production together with (|j) forces the gluino mass to be too large. It should be 
noted that M3 as given by (Q) is less than the physical pole mass of the gluino |2g] by an amount 
which is often quite substantial, especially if the squarks are heavier. This effect makes it even 
less likely that processes involving gluino production can compete with the chargino/neutralino 
processes already discussed. The rate for associated production of Cj or Ni with a gluino 
or squark overtakes the (gg + gq + qq) production rate when trig exceeds roughly 400 GeV, 



depending on the squark masses [29]. In this regime, however, even the sum of all processes 
involving gluino or squark production should be small compared to those from chargino and 
neutralino production, unless rriq is significantly less than m g . It therefore seems unlikely that 
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gluinos or squarks can be involved in the discovery process. 

E. Comments on the interpretation of the CDF ecyy^T event 



In this section we have emphasized the power of the Tevatron in setting exclusion limits in 
the gravitino LSP scenario. Of course, in this context we must mention that at least one event 



[12] of this general type has been observed at CDF. This event has an energetic electron and 
positron, two energetic photons each with \n\ < 1, and large (> 50 GeV) ]fcp. It easily passes 
the cuts defining our signal above. 

The most obvious candidate process for this event is selectron pair production. As has 
already been discussed in 14], and recently in some more detail in 0], one can attempt to 



explain the event either in terms of e R e R pair production or e L e L pair production. From the 
kinematic information, we found in [14] that in either of these two cases, one has rough bounds 



nig > 80 GeV and 38 < < 100 GeV. If one assumes gaugino mass unification, the lower 
bound obtained here for is far weaker than the lower bound established above from non- 
observation of chargino and neutralino events at the Tevatron. As has been emphasized recently 



in [11], the energetic electrons in the event seem to indicate a significant mass difference between 
rrifj and mg, in order to have sufficiently energetic electrons with a high enough probability to 
explain the event. 

Right-handed selectrons have a lower production cross-section than do left-handed selec- 
trons for a given mass, as can be seen from Fig. and this seems to perhaps favor the idea 
that the pair-produced selectron was left-handed. However, since we are forced to calculate 
the probability of this single event using Poisson statistics, this argument is not on very solid 
footing. For example, if the three right-handed sleptons are degenerate in mass, then the cross- 
section to produce any pair £ R £ R is of course 3 times larger than the rate for the pair e R e R which 
could explain the event actually seen. It is not completely clear which of these rates should 
be taken in assessing the likelihood of a single observed event. In any case, the rate before 
cuts for £ R £ R production in 100 pb _1 is about 1 event for m« = 125 GeV and 1/2 event for 



= 145 GeV. The rates after cuts are significantly less because of acceptances (see Fig. 10), 
but it seems possible that the observed event is an upward fluctuation in the e R e R production 
process, even with nig R — sufficiently large to explain the observed kinematics. 

Conversely, although the rates for left-handed selectron production are larger, one must 
also note that in this interpretation the expected number of events with two leptons is always 
considerably less than for one or zero leptons (see Fig. [ll]). The reason for this is that the 
production cross-section from i>v and vi L is necessarily larger than for £ L £ L because of the 
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kinematics dictated by the sum rule (21). It might therefore be viewed as problematic that 
the single observed event has two leptons. The limited acceptance for leptons only exacerbates 
this problem. Nevertheless, it again seems to be not entirely out of the question that the event 
could be due to e L e L production. 

There is another very interesting possibility, however, illustrated by model 4 above, that 
the event could be due to chargino pair-production. We note that depending on parameters, 
the chargino pair-production cross section remains sufficiently large to give 1 event (after cuts) 
in 100 pb _1 up to at least m,Q = 200 GeV. Now, each produced chargino can decay into either 
IvyG or qq^G. If the decay is dominantly through VK-boson or squark exchange, then one 
might expect the ee^lpT event to be accompanied by many more events with jets or significant 
hadronic activity in addition to However, this can be avoided if C\ has kinematically 

allowed 2-body decays into left-handed sleptons and not into WN±: 



> mo > mfa ; (32) 
^ < raw- (33) 



v(Iq — m 

In this case, the chargino should decay as either C\ — > vl — » lvN\ — > Ij^t or (if kinematically 
allowed) C\ — > t h v — > lvN\ — > £j$t- Both of these decays have the same signal (with different 
kinematics), so that the signal for chargino pair production will be £ + i'^jj^T in this case, 
with a nearly 100% branching fraction before cuts. Even though sneutrinos are lighter than 
charginos, the chargino production cross-section can be much larger. Of course, it should be 
noted that the rate for different flavor leptons in this case is twice that for like-flavor leptons. 
Still, we find that it is possible to obtain kinematics and rates after cuts which could explain 
the CDF event. The kinematics of the event together with the cross-section for chargino pair- 
production evidently favor mg > 125 GeV in this case. 

We should also mention that N1N2 production can give an eeyy^T signal, but it is very 
difficult to reconcile this possibility with the observed event, because of the large invariant mass 



of the ee pair in the event [12|. 



4. Supersymmetry with a light gravitino at LEP 



The LEP collider at CERN will probe some of the parameter space which is not yet 
excludable by the current 100 pb _1 data sample collected at the Tevatron. However, it is 
important to take into account the results of Section 3 when assessing the discovery potential 
of the various LEP upgrades. At least within the context of gaugino mass unification, we 
have found that the lightest neutralino mass can be bounded from below by 70 GeV, based on 
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the exclusion capability of the current 100 pb data. Similarly, the lighter chargino mass is 
bounded below by 125 GeV, and even when (Q) is not assumed, one has > 100 GeV for 
mft > 50 GeV. Therefore, it is immediately clear that one cannot hope to observe chargino 
pair production at any of the LEP upgrades considered here in the gravitino LSP scenario 
with our assumption that N\ — > G7 always occurs within the detector. Furthermore, the 
second lightest neutralino should also not be kinematically accessible at LEP even in N1N2 
production, at least in the case that gaugino mass unification (|j) holds. The reason for this is 
that m ( j i > 125 (as required by the Tevatron data) and < 95 (as required for accessibility 
in e + e~ collisions with ^/s = 190 GeV) forces one into a region of parameter space with rather 
large \fj,\ and gaugino-like N± and N%, so that rrifj + > 210 GeV. Therefore, it is clear 
that in the chargino/neutralino sector, LEP190 can only hope to observe N\N\ production 
with signature 77^. Likewise, the existing Tevatron data makes it impossible for a light stop 
(or other squark) to be accessible at LEP with our assumptions. There is a still a possibility 
to observe slepton pair production since, taking into account efficiencies, the Tevatron cannot 
set exclusion limits on slepton masses which are significantly stronger than the indirect one 
following from > > 70 GeV. Therefore there is a narrow range of N\ and slepton 

masses from no less than 70 GeV up to less than 95 GeV which can be probed at LEP with 
y/s < 190 GeV. 

We begin by considering N\N\ production in e + e~ collisions, which leads to events with 
two acoplanar photons and large missing energy. (A similar study for the NLC has recently been 
made ||].) The energy distribution of photons produced^ in such events is flat, with endpoints 

-EWn < -E71 j -E72 < -Emax (34) 
^max.min = j{Vs ± \J S ~^ m \)- ( 35 ) 

The two photon energies in each event vary over this range independently, providing a very sim- 
ple characteristic kinematic signature. The missing energy in each event is bounded according 
to 2E m { n < ]fi < 2£^ max and is peaked at -Ebeam = y/s/2. Two further corollaries are that the 
distribution of E^ + E^ 2 is the same as that of fi, and that the energy distribution of the more 
(less) energetic photon observed in each event rises (falls) linearly with energy. The numerical 
bounds on photon energies in the N\N\ signal are, for the various LEP upgrades: 

20 GeV < £ 71 ,£ 72 < 60 GeV (y/s = 160 GeV) (36) 
18 GeV < E yi ,E 72 < 70 GeV (y/s = 175 GeV) (37) 
16 GeV < E yi ,E J2 < 80 GeV (y/s = 190 GeV) (38) 

for a lower bound =70 GeV. For masses nearer threshold, the range of photon energies of 
course becomes narrower around \/s/A in each case. Thus the lower bound on Ni mass from 
3 We neglect final state interference effects throughout the following discussion. 
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the Tevatron ensures that the NiNx signal at LEP will automatically pass appropriate cuts on 
soft photons. This will be useful below in our discussion of cuts and backgrounds. 

Several factors affect the production cross-section for NiNi at LEP. Since in the accessible 
parameter space N\ is essentially forced to have a large gaugino component, the s-channel Z 
boson exchange contribution is suppressed. If sleptons are light, the diagrams with slepton 
exchange will dominate. The diagrams with e R exchange are usually far more important, 
because the ee R N\ coupling is larger than the ee L N\ coupling. The N\N\ production cross- 
section is quite sensitive to the selectron masses, even if the selectrons themselves are not 
accessible at LEP. As a result, the discovery reach is always within a few GeV of the kinematic 
limit, but for no value of Trijj can one clearly guarantee discovery at any of the LEP upgrades, 
because of low observable rates for large mg R . At LEP160, the cross-section is always less than 
0.2 pb for models (with gaugino mass unification) not excludable at the Tevatron, and is less 
than 0.1 pb for > 75 GeV. These are optimistic upper bounds, and the cross-sections for 
less favorable parameters can be much smaller. This leaves open only the possibility of perhaps 
a few events at LEP160 with 25 pb -1 per experiment, for a narrow mass range, optimistically 
70 GeV < < 77 GeV. As we will remark below, there is also a non-trivial background for 
such events, so that an unambiguous discovery will require a certain amount of luck. 

The prospects for discovery (or confirmation) are clearly much brighter at LEP 190 with 500 
pb -1 per experiment, both because of the kinematic reach and the greater luminosity. In Figure 



15, we show a scatter plot of the total N±N\ cross-section at y/s = 190 GeV. Each point on this 
plot corresponds to a set of model parameters which plausibly could have avoided detection 
at the Tevatron with the current integrated luminosity, based on the results of the previous 
section. To illustrate the dominance of the e R exchange diagrams, models with mg R < 175 GeV 
are denoted by crosses, while those with 175 GeV < mg H < 500 GeV are denoted by dots. (The 
CDF eejjlpT event could perhaps be explained by e R e R production in models of the former 
category.) The gaugino mass unification condition (Q) is assumed here. We have taken into 
account initial state radiation effects which imply a small ( < 10%) reduction in the signal; 
it should be noted that such effects are larger when Z boson exchange dominates because of 



radiative return to the Z peak. To these cross-sections one must apply detector cuts, e.g. |30| : 

|cos6> 7 | < 0.95, (39) 
(p T )j > 0.065 £ beam , (40) 

for each photon. The discovery reach with 500 pb -1 extends up to within a few GeV of the 
kinematic limit. Clearly the presence of a light e R provides much more favorable discovery 
prospects. However, there is no guarantee of discovery of a light N\ if mg R is large, even with 
this amount of luminosity, and even for the most favorable kinematics. 
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Figure 15: Production cross-section for N±Ni at LEP with -y/s = 190 GeV. Each point represents 
a set of model parameters which may not be excludable at the Tevatron with 100 pb _1 . The 
crosses (dots) represent models with mg R less (greater) than 175 GeV. 

We now turn to the question of backgrounds. The ordinary QED process e + e~ — > 77 
production has a large cross section, but can easily be discriminated against with a cut on 
missing energy or equivalently E~ < 0.8E^ caum for each photon. The most important physics 
backgrounds for the 77^ signal come from 77z/j!7j (i = e, fi, r) with two separately gauge- 
invariant sets of diagrams: 

A) e + e~ — ► 77ZW with Z — > vjpi. (3 Feynman diagrams); 

B) e + e~ — > 77f e I7 e through virtual W^-boson exchange. (7 Feynman diagrams). 

We have computed these backgrounds using CompHEP |3l]] , a specialized package for au- 
tomated calculation of high-energy elementary particle processes, with results fed into BASES, 
a Monte Carlo phase-space integration program. The processes e + e~ — > 77^^ and e + e _ — * 
77z/ T F r receive contributions only from the type A diagrams. At = (160, 175, 190), they each 
contribute (49,37,30) fb to the background for 77^ after the cuts ©, @). When the final 
state is 771^ I7 e , one must take into account a significant interference between the diagrams of 
types A and B. The diagrams of type A clearly dominate in the kinematic regime characterized 
by a missing invariant mass Mi nv ; s very close to Mz- [Here Mj 2 nvis = (p e + + p e - — p 7l — p 72 ) 2 -] 
In that regime, the interference with the type B diagrams are a small perturbation and in any 
case only affect 1/3 of the background. For slightly larger Minvisj however, the type B diagrams 
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do have a substantial interference with the off-peak type A diagrams. The overall effect is one 
of constructive interference, but the sign is not definite for all kinematic configurations. At 
y/s = (160, 175, 190), we find (61,49,42) fb for e+e" -> ^u e V e after the cuts ©, ©. 

Since the 77Z/F backgrounds have larger support for relatively soft photon energies, one 
can reduce them somewhat by imposing the cut 

0.2 < E 7 /E beam < 0.8 (41) 

on each photon; the upper limit easily eliminates the e + e~ — > 77 process as we have already 
mentioned. The cut (^) has little or no effect on the signal, as can be seen from (p4])-(|38|). After 
imposing this cut in addition to the detector cuts (|39|), fl4C|), we find a remaining background 
at y/s = (160, 175, 190) of (29, 24, 21) fb from jju e u e , and (27, 22, 18) fb from each of 771^ 
and 77^ r I7 T . 

In order to more strongly reduce the backgrounds we can impose a cut on the missing 
invariant mass of 

10 GeV < Mi nvis < 80 GeV . (42) 

The upper limit is to avoid the 77Z/F physics backgrounds, while the lower limit eliminates a 
potentially large (several hundred fb) detector background following from e + e~ — > 77(7) with 
one photon lost in the beam direction or in an insensitive part of the detector. (The part of this 
background due to photons lost in the beam direction is also substantially reduced by imposing 
a lower bound cut on fa-) The signal vanishes at the endpoints of the distribution Mi nv ; s = 0, 
2-Emax and is broadly distributed in between. The greater part of the signal will always pass all 
of the cuts, although a significant part of the signal will necessarily have to be eliminated by 
the cut ([42D . After imposing this cut in addition to (p9|), (|40|), the total 771/F backgrounds at 
y/s = (160, 175, 190) GeV are only (1.8, 1.3, 1.0) fb respectively. Finally, imposing the cut @ 
on top of these cuts reduces the background to a completely negligible level. 

The distribution in Mi nv i s for signals and backgrounds at y/s = 190 GeV are shown in 



Figure 16. In this figure we have arbitrarily chosen a total signal (before cuts) of 50 fb, with 
Ni masses of 75 and 90 GeV. The differential cross-sections shown are after the detector cuts 
(|3"9D, (f40|) and photon energy cuts (gl]). The total 77^z/ background shown amounts to 56 fb, 
but is reduced to a negligible level by the Mi nv i s cut. Note, however, the significant overlap in 
invariant missing mass for the backgrounds and the signals. The signal distribution in Mi nv ; s is 
broadly peaked below the 80 GeV cut, and vanishes near Mi nv ; s = 0. We conclude that even in 
the worst-case kinematic situation, the efficiency for detecting N±Ni should exceed 50% after 
cuts at LEP190. Thus a 40 fb signal before cuts should provide a 10 event discovery after 
cuts with 500 pb _1 . By comparing with Fig. 15, we conclude that LEP190 should be able to 
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Figure 16: Distribution of the missing invariant mass in 77^ events at LEP with i/s = 190 
GeV. Angular and photon energy cuts have been applied as described in the text. The lighter 
solid line is the remaining total background (56 fb) for all three neutrino species. The signals 
for mft = 75 and 90 GeV are the solid and dashed lines, respectively, with an arbitrary choice 
of 50 fb for the signal before cuts in each case. 

unambiguously observe NiN\ production for up to at least 85 GeV if mg R < 175 GeV, 
assuming gaugino mass unification. The exclusion capability decreases for larger m e - fl , however. 
The 771/77 background is more problematic at y/s = 160 GeV with 25 pb~ x or yfs = 175 GeV 
with 10 pb" 1 , where only at most a few signal events are expected, and the signal distribution 
in Mi nv i s again overlaps with the Z-boson peak. 



Conversely, Figure 17 shows the distributions for photon energies at yfs = 190 GeV for 
the backgrounds, before and after the Mi nv i s cut. All of the 77Z/F contributions have been 
included here. The two distributions correspond to the more and less energetic photon in each 
background event, after the detector cuts (3S), (40). After imposing in addition the cuts ( |4TD 
and ([42]), the total 77Z/F background is reduced to a fraction of a femtobarn. The signal from 
N\N\ production (not shown) is characterized by a linearly rising (falling) distribution for the 
more (less) energetic photon in each event, with endpoints -Emm and E max as given above. (Note 



the log scale in Fig. 17. 



We now turn to the question of slepton pair production signals at LEP. In general, slep- 
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Figure 17: Distribution of photon energies for 771/F backgrounds at LEP with yfs = 190 
GeV. Detector cuts have been applied as described in the text. The solid (dashed) line is the 
distribution for the more (less) energetic photon 72 (71) in each event. The dotted (dot-dash) 
lines are the same distributions after the cut on Mi nv ; s described in the text. 

ton masses up to within a few GeV of the kinematic limit should lead to visible signals with 
500 pb _1 at LEP190. If right-handed sleptons are kinematically accessible, one finds that the 
cross-section for e R e R production is generally somewhat larger than those for each of fi R jl R and 
f R f R , because of the positive contribution of diagrams with t-channel exchange of gaugino-like 
N\. The pair production of electron sneutrinos can be very strongly suppressed, because of 
destructive interference from chargino exchange, even with C\ required to be heavier than 125 
GeV. Fortunately, pair production of muon and tau sneutrinos does not suffer this suppression, 
and those cross-sections are always large up to within a GeV or two of the kinematic limit. 
Because of the sum rule (^H), it seems quite unlikely that pair production of left-handed se- 
lectrons can be a discovery process at LEP in the gravitino LSP scenario considered in this 
paper. However, the cross-section for e L e R can be even larger than for e R e R production when 
both are kinematically accessible, because of a large contribution from exchange of gaugino-like 
neutralinos. 

It is important to note that for the (quite narrow) range of masses which are accessible at 
LEP and which cannot already be ruled out at the Tevatron, each slepton has one (and only 
one) allowed 2-body decay mode, namely I — ► N±£. This decay is never strongly suppressed 
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because Ni always has a significant gaugino component. Therefore, charged slepton production 
will essentially always give rise to the signal while sneutrino production, like N\Nx 

production, can give rise only to 77^. The leptons appearing in jj£ + l~]fi events from charged 
slepton pair production at LEP should necessarily be quite soft, because there cannot be a 
large mass difference between the slepton and N\ and the sleptons cannot have a large boost. 
However, the SM backgrounds for such processes are extremely small. Taking into account the 
cuts (|39|), (|40|), one finds that 77ZZ production is always below threshold at LEP160, LEP175, 
and LEP190, while j^yWW is only above threshold at LEP190. Using CompHEP we have found 
that the latter process only contributes about 0.1 fb to the ^77^ background at yfs = 190 
GeV. There is also a background for ee^f and finjjlfi from the process e + e~ 77 with 

Z — > t + t~ and leptonic r decays, but this is very small. Similarly, the photons produced in 
sneutrino pair production should be softer than those found in N\N\ events. Since slepton 
interactions are not expected to exhibit significant flavor violation, we can conclude by noting 
that the signatures for the gravitino LSP scenario at LEP are always 77^ (from N\N\ and vv 
production) and e + e~77^, fi + fi~j^/]fi, and r + r~77|5 (from charged slepton production). 



5. Concluding Remarks 



In this paper, we have studied discovery signals for supersymmetry with a gravitino LSP 
at the Tevatron and at LEP. If the decay N\ — * 7G occurs within the detector, then super- 
symmetric phenomenology at colliders will have a very bright future. Indeed, the existing 
Tevatron data of 100 pb _1 should allow the exclusions > 70 GeV and mg > 125 GeV 
at least in models obeying the gaugino mass unification condition @). For > 50 GeV, it 
should be possible to exclude mg < 100 GeV in a model-independent way. The reach is much 
higher. These results rely on the fact that every supersymmetric event contains two potentially 
detectable energetic photons and $t, yielding a high (up to 30%) detection efficiency. We em- 
phasize that this efficiency is not expected to be significantly reduced by small mass splittings 
between charginos and neutralinos, since both the photon energies and the $t depend only on 
the mass and boost of the N\. If the single CDF ee^^IpT event is an example of such an event, 
then it is not unlikely that an upgraded Tevatron with yfs = 2 TeV and > 2 fb -1 of data can 
establish a discovery. In any case, Tevatron upgrades will continue to make strong inroads into 
the parameter space of the gravitino LSP scenario. The reach and exclusion capability can be 
estimated for Tevatron upgrades using Figs. || and |2| and assuming efficiencies > 15% as found 
in Figs. H and ||. (Note that the increase in yfs from 1.8 TeV to 2 TeV makes such estimates 
conservative.) 
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It is quite possible that LEP190 with 500 pb _1 of data can make the discovery if the lightest 
neutralino is kinematically accessible. At least in the case of models not already excludable by 
the Tevatron with gaugino mass unification, only N±Ni and slepton pair-production can be ex- 
plored at LEP, with possible signals 77-^ and t r t~^'$. We found that appropriate cuts on the 
missing invariant mass and on photon energies can reduce the 77^ backgrounds to a negligible 
level while keeping intact at least 50% of the signal even in the worst kinematic situation. The 
discovery reach extends to within a few GeV of the kinematic limit. An important factor in the 
e + e _ — > N\Ni —* 77^ search is the mass of the right-handed selectron. If rrie R ^ 175 GeV, 
then there should be at least 10 events after cuts in 500 pb~* at y/s = 190 GeV for ^85 
GeV, but the rate can be much lower for larger mg R - In this sense, any exclusion limits will 
be dependent on assumed upper bounds for mg R . If e R is light, then LEP160 and LEP175 can 
observe a few events. 

Although we have not studied future colliders here, it seems clear that both the Large 
Hadron Collider and Next Linear Collider will be very effective discovery machines if the de- 
tectors have good efficiency for detecting isolated energetic photons and fir- If the reported 
ee77-^T CDF event is interpreted as slepton or chargino production, it seems essentially certain 
that the NLC will detect supersymmetric events, and that the LHC also will if the detectors 
are sufficiently good with photons and $t- 

In general, the ability of the present Tevatron data sample to bound the gravitino LSP 
scenario emphasizes the importance of photon detection. This component should not be ignored 
in future detector design. Also, it would be useful to have photon pointing information, in case 
the Ni — ► 7G decay length is macroscopic. As we mentioned in section 2, it is certainly 
possible that this is so, leading to more single photon events than diphoton events. In that 
case, one can imagine discovering super symmetry using the usual well-known discovery signals 
of the neutralino LSP scenario, supplemented by a fraction of these events with one additional 
energetic photon. Measuring this fraction would provide a powerful piece of information in 
disentangling the signal, especially if it can be combined with measurements of the Ni decay 
length. 

The gravitino LSP possibility also provides a rich area for theoretical explorations. This 
scenario necessarily implies a low scale of spontaneous supersymmetry breaking, but the precise 
mechanism for this breaking and its communication to the fields of the MSSM remain open 
questions ||. The connections between supersymmetry breaking and the ^ parameter and other 



aspects of electroweak symmetry breaking are also interesting questions [32|. We should also 
mention that if the gravitino is the LSP, then the lightest neutralino is of course no longer a cold 
dark matter candidate. It remains to be seen if one can obtain a viable dark matter scenario; 
for a recent proposal see p3[. It seems clear that such issues merit further investigation. 



34 



Acknowledgments 



We are indebted to Bob Blair, Marty Einhorn, Henry Frisch, Tony Gherghetta, Joey Hus- 
ton, and Eve Kovacs for helpful discussions. This work was supported in part by the U.S. 
Department of Energy. 



References 

[1] P. Fayet, Phys. Lett. B 70 (1977) 461; Phys. Lett. B 86 (1979) 272; Phys. Lett. B 175 
(1986) 471; and in "Unification of the fundamental particle interactions", eds. S. Ferrara, 
J. Ellis, P. van Nieuwenhuizen (Plenum, New York, 1980) p. 587. 

[2] M. Dine, W. Fischler, M.Srednicki, Nucl. Phys. B189 (1981) 575; S. Dimopoulos, S. Raby, 
Nucl. Phys. B192 (1981) 353; M. Dine, W. Fischler, Phys. Lett. B 110 (1982) 227; M. Dine, 
M. Srednicki, Nucl. Phys. B202 (1982) 238; L. Alvarez-Gaume, M. Claudson, M. Wise, 
Nucl. Phys. B207 (1982) 96; C. Nappi, B. Ovrut, Phys. Lett. B 113 (1982) 175; M. Dine, 
W. Fischler, Nucl. Phys. B204 (1982) 346; S. Dimopoulos, S. Raby, Nucl. Phys. B219 
(1983) 479. 

[3] For recent examples see M. Dine, A. E. Nelson, Phys. Rev. D48 (1993) 1277; M. Dine, 
A. E. Nelson, Y. Shirman, Phys. Rev. D51 (1995) 1362; M. Dine, A. E. Nelson, Y. Nir, 
Y. Shirman, Phys. Rev. D53 (1996) 2658. 

[4] A. E. Nelson, N. Seiberg, Nucl. Phys. B416 (1994) 46; J. Bagger, E. Poppitz, L. Randall, 
Nucl. Phys. B426 (1994) 3. 

[5] T. Bhattacharya, P. Roy, Phys. Lett. B 206 (1988) 655; Nucl. Phys. B328 (1989) 481. 

[6] H. Pagels, J. R. Primack, Phys. Rev. Lett. 48 (1982) 223; T. Moroi, H. Murayama, M. Ya- 
maguchi, Phys. Lett. B 303 (1993) 289. 



[7] For a review, see H. Baer et al., "Low-energy supersymmetry phenomenology", hep- 



ph/9503479 



[8] N. Cabibbo, G. R. Farrar, and L. Maiani, Phys. Lett. B 105 (1981) 155; M. K. Gaillard, 
L. Hall, I. Hinchliffe, Phys. Lett. B 116 (1982) 279; J. Ellis, J. S. Hagelin, Phys. Lett. B 
122 (1983) 303; D. A. Dicus, S. Nandi, and J. Woodside, Phys. Lett. B 258 (1991) 231. 



[9] D. R. Stump, M. Wiest, C. P. Yuan, preprint |hep-ph/9601362 



[10] K. Kiers, J. N. Ng, G. Wu, preprint |hep-ph/ 96043381 



[11] S. Dimopoulos, S. Thomas, J. D. Wells, preprint hep-ph /9604452 



[12] S. Park, "Search for New Phenomena in CDF", 10 th Topical Workshop on Proton-Anti- 
proton Collider Physics, edited by Rajendran Raja and John Yoh, AIP Press, 1995. 



35 



S. Dimopoulos, M. Dine, S. Raby, S. Thomas, Phys. Rev. Lett. 76 (1996) 3494. 

S. Ambrosanio, G. L. Kane, G. D. Kribs, S. P. Martin, S. Mrenna, Phys. Rev. Lett. 76 
(1996) 3498. 

E. Cremmer, S. Ferrara, L. Girardello, A. van Proyen, Nucl. Phys. B212 (1983) 413. 

H. E. Haber, G. L. Kane, Phys. Rep. 117 (85) 75; J. F. Gunion, H. E. Haber, Nucl. Phys. 
B272 (1986) 1; Erratum ibid B402 (1993) 567. 

J. M. Cornwall, D. N. Levin, G. Tiktopoulos, Phys. Rev. D10 (1974) 1145, Erratum ibid 
Dll (1975) 972. 

D. A. Dicus, S. Nandi, J. Woodside, Phys. Rev. D41 (1990) 2347. 

M. Drees, J. Woodside in Super symmetry, Proc. of ECFA Workshop on the Large Hadron 
Collider, 1990. 

M. Drees, J. Ellis, P. Jetzer, D. W. Sciama, Phys. Lett. B 220 (1989) 586 

H.U. Bengtsson, T. Sjdstrand, Comp. Phys. Comm. 46 (1987) 43; T. Sjostrand, Comp. 
Phys. Comm. 82 (1994) 74; S. Mrenna, "Simulating Supersymmetry With PYTHIA 5.7 
and JETSET 7.4," CITHE-68-1987, July 1995. 

See for example H. Baer, C.-h. Chen, C. Kao, X. Tata, Phys. Rev. D52 (1995) 1565; 
S. Mrenna, G.L. Kane, G.D. Kribs, J.D. Wells, Phys. Rev. D53 (1996) 1168. 

We thank Bob Blair (Argonne) for providing information on the CDF diphoton cross- 
section measurements. 

We thank T. Sjostrand (Lund) for this suggestion. 

H. Baer, C.-h. Chen, F. Paige, X. Tata, Phys. Rev. D49 (1994) 3283. 

A. Djouadi et al. Nucl. Phys. B349 (1991) 48; M. Boulware, D. Finnell, Phys. Rev. D44 
(1991) 2054; J. D. Wells, C. Kolda, G. L. Kane, Phys. Lett. B 338 (1994) 219; D. Garcia, 
J. Sola, Phys. Lett. B 354 (1995) 335; G. L. Kane, R. G. Stuart, J. D. Wells, Phys. Lett. B 
354 (1995) 350; A. Dabelstein, W. Hollik, W. Mosle, [hep-ph/9506~25l| ; P. H. Chankowski, 
S. Pokorski, Phys. Lett. B 366 (1996) 188; J. Ellis, J. Lopez, D. Nanopoulos, Phys. Lett. 
B 372 (1996) 95; J. D. Wells, G. L. Kane, Phys. Rev. Lett. 76 (1996) 869; E. H. Simmons, 
Y. Su, |hep-ph/9602267| ; P. H. Chankowski, S. Pokorski, |hep-ph/96033l0| . 

[27] D.R. Claes, 10 th Topical Workshop on Proton-Anti-proton Collider Physics, edited by 
Rajendran Raja and John Yoh, AIP Press, 1995, FERMILAB-CONF-95/186-E. 

[28] S. P. Martin, M. T. Vaughn, Phys. Lett. B 318 (1993) 331. The corresponding (smaller) 
corrections for charginos and neutralinos are given in D. Pierce and A. Papadopoulos, 
Phys. Rev. D50 (1994) 565, Nucl. Phys. B430 (1994) 278. 

[29] H. Baer, C.-h. Chen, C. Kao, and X. Tata, Phys. Rev. D52 (1995) 1565. 



36 



[30] F. Boudjema, B. Mele (Conveners), E. Accomando et al. [The Standard Model Process 
Group], Proc. of the Workshop "Physics at LEP2", G. Altarelli, T. Sjdstrand, F. Zwirner 



(eds.), CERN-Report 96-01 (1996), [hep-ph/9601224 . 



[31] E. E. Boos, M. N. Dubinin, V. A. Ilin, A. E. Pukhov, V. I. Savrin, "CompHEP: Spe- 
cialized package for automatic calculation of elementary particle decays and collisions", 



hep-ph/9503280 , and references therein. 



[32] G. Dvali, G. F. Giudice, A. Pomarol, |hep-ph/9603238| . 



[33] S. Borgani, A. Masiero, M. Yamaguchi, hep-ph/9605222 . It is not clear that the require- 
ments of this paper can be satisfied if the CDF eerpyjfir event is interpreted as superpartner 
production. 



37 



